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We perform a thorough theoretical analysis of resonance effects when an arbitrarily polarized plane 
monochromatic wave is incident onto a double periodically modulated metal film sandwiched by two 
different transparent media. The proposed theory offers a generalization of the theory that had been 
build in our recent papers for the simplest case of one-dimensional structures to two-dimensional 
ones. A special emphasis is placed on the films with the modulation caused by cylindrical inclusions, 
hence, the results obtained are applicable to the films used in the experiments. We discuss a spec- 
tral composition of modulated films and highlight the principal role of "resonance" and "coupling" 
modulation harmonics. All the originating multiple resonances are examined in detail. The trans- 
formation coefficients corresponding to different diffraction orders are investigated in the vicinity 
of each resonance. We make a comparison between our theory and recent experiments concerning 
enhanced light transmittance (ELT) and show the ways of increasing the efficiency of these phe- 
nomena. In the appendix we demonstrate a close analogy between ELT effect and peculiarities of a 
forced motion of two coupled classical oscillators. 

PACS numbers: 42.25.-p, 73.20.Mf, 78.20.Ci, 78.20.Bh 



I. INTRODUCTION 

It is exactly tangible advances in structuring metals 
on nanoscales that give rise to a grate amount of exper- 
imental and theoretical works in the field of plasmonics, 
which examine resonance optic effects caused by surface 
plasmon polaritons (SPPs) excitation in structured con- 
ducting or semiconducting media 1-3 . These quasi two- 
dimensional (2D) electrodynamic objects have been thor- 
oughly studied in the solid state physics, physics of sur- 
faces and diffraction optics over the late few decades 4 . An 
SPPs are electromagnetic surface waves, coupled to the 
collective electron excitation, they attract a great deal of 
attention due to their unique possibility of light localiza- 
tion and essential enhancement of the electric field near 
the surface. 

A considerable interest in SPPs stems from the lat- 
est experiments on ELT phenomena. Since 1998, after 
observation by Ebbesen et.al. 5 ' 6 the violation of Bethe's 
approach 7 to the diffraction by subwavelength periodic 
hole arrays in real metal films, this has been the sub- 
ject of many studies. Until the recent moment one of 
the most generally recognized explanations of the ELT 
through subwavelength periodic hole arrays is the exci- 
tation of SPPs. Most of the authors (especially theoreti- 
cians), who hold on to the SPP conception of ELT, as- 
sume that the field enhancement results from "interface" 
SPPs, being either single-boundary localized (which is a 
common case for a nonsymmctric dielectric surrounding 
of the film) or double-boundary localized for the sym- 
metric surrounding. A periodic hole array plays a role of 
a coupler between the SPPs existing at the boundaries 



of the conducting film and the incident light. In this 
context the crucial point is the surface periodicity. The 
periodicity caused by other factors, such as corrugation 
or periodic modulation of the medium electromagnetic 
properties, etc., also result in the light-matter interaction 
resonance features, in particular, the ELT effect. Note 
that comparing the ELT transmittance peaks positions 
to those caused by SPP excitation in different diffraction 
orders undoubtedly points out a significant SPP role, see 
numerous experiments, e.g. Refs. 5,6,8-14. 

It should be noted that the observed ELT effects 
are strongly dependent upon the film surrounding. As 
the film having the subwavelength hole array is sur- 
rounded by the dielectrics with the same dielectric con- 
stants (for instance, a free-standing film) the excited 
SPPs are double-boundary- localized, that is, the field is 
enhanced at both faces of the film, see experiments 8 ' 9 . 
In this geometry the ELT in zeroth diffraction order is 
much more pronounced as compared to the nonsymmet- 
rical film surrounding, when the excited SPPs are single- 
boundary-localized, that is the field can be enhanced at 
a single face of the film. Generaly, this corresponds to 
the films deposited onto the quartz substrate. Both of 
these cases have been so far studied by using numerical 
methods 8,9,14-19 . It was also shown numerically and ex- 
perimentally that the ELT occurs for periodical conduct- 
ing structures without holes. Basically this is quite evi- 
dent, since the type of periodicity does not play a crucial 
role for the excitation of the interface SPPs. These were 
chiefly the structures with relief corrugations of the film 
faces both for symmetrical 20 and nonsymetrical 21-23 di- 
electric surrounding and the structures with periodically- 
located dielectric pillars 24 . 
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Several authors have developed an analytical approach 
which qualitatively describes the ELT. They have stud- 
ied the diffraction by the film with one-dimensional (ID) 
periodically modulated dielectric permittivity 25-30 . In 
these works a study has been made in the simplest case 
of a strictly normal incidence onto the symmetrically sur- 
rounded film with harmonic modulation of the film per- 
mittivity. The exception are the papers 27,30 , where the 
theory is generalized to the nonsymmetric surrounding. 
The authors of these works have described the zeroth- 
order transmittance dependence upon the parameters in 
case of SPP excitation in the first diffraction order. In 
contrast to them, in a certain sense, we have provided a 
more general analytical insight 31-33 , by solving the prob- 
lem of vector diffraction by the film with nonsymmetric 
and symmetric dielectric surrounding for arbitrary mod- 
ulation Fourier spectra of the permittivity of the metal at 
an arbitrary incident angle and at an arbitrary polariza- 
tion of the incident light. The advantage of our analytical 
treatment is that we have described not only the first- 
order resonances, but have given a classification of the 
resonances corresponding to single or double-boundary- 
localized SPPs excitation in single or multiple diffraction 
orders. Also, we have examined a nonzcroth-order ELT 
that was observed in the experiments Refs. 22,34. Men- 
tion also paper 35 , where the alternative analytical ap- 
proach was suggested to describe the light transmittance 
through metallic nanoslit structures. 

Yet another important class of plasmonic structures 
possessing interesting diffraction resonance phenomena 
are metallic nanoparticle arrays 12,36-38 . As known, they 
can support the so-called localized SPP resonances which 
arc strongly dependent upon the shape of the individ- 
ual particles 38 . It is believed that the shape resonances 
similar to the localized ones, may affect the ELT in sub- 
wavelength hole arrays. The influence of the hole or/and 
nanoparticle shape was studied experimentally in 39-41 . 
The study of the conducting films containing periodic lat- 
tices of dielectric nanoparticles and voids have been made 
in Refs. 42-44. However, we are not aware of the study 
on the experimental or theoretical optical properties of 
metallic nanoparticle arrays immersed into a conducting 
film. Similarly to a nanohole array, such a nanoparticle 
array should display ELT with the wavelength spectra 
depending strongly upon the nanoparticle shape. 

In the present paper we go on investigating the reso- 
nance optical effects by generalizing the earlier developed 
analytical treatment 31,32 of the conducting films with 2D 
modulation. We consider the vector diffraction problem 
for periodically located metallic inclusions in the metal 
film with an arbitrary dielectric surrounding. Since the 
inclusions are presumed to be entirely imbedded into the 
film (the faces of the film being flat), the resonances of 
our system correspond to the excitation of purely inter- 
face SPPs, the localizcd-SPP resonances do not exist. 
Along with the fact that such structures arc of inter- 
est by themselves, they may be considered as a model 
of subwavelength hole arrays; and the approach devel- 



oped is appropriate for other periodical structures (say, 
for corrugated films). One more substantial aspect of our 
approach is that we can easily describe the polarization 
of light transmitted or/and reflected by the 2D periodical 
structure. The polarization properties of subwavelength 
hole arrays is also a matter of interest, and they have 
been intensively studied in recent years 33,45-49 . 



The paper is arranged as follows. Following the In- 
troduction, in Section II, we describe a general approach 
to the problem of resonance light diffraction by a 2D 
periodically-modulated conducting film for the conical 
mount case and for an arbitrary polarization of the in- 
cident light. We accentuate the fact that the form of 
the periodically-located inclusions have an impact upon 
the Fourier spectra of the periodical structure and this, 
in its turn, influences the excited SPPs considerably. In 
Section III we examine the excitation of single-boundary- 
localizcd SPPs in the nonsymmetrically-sandwichcd film, 
considering both the entire transmittance and reflectance 
spectra and a more detailed analysis of different reso- 
nances in their close vicinity. We give an explanation of 
the recently observed polarization dependence upon the 
hole shape 39 from the view point of Fourier spectra of 
the periodical structure. Besides, we compare our cal- 
culations with other recent experiments. Section IV is 
devoted to the resonance effects caused by the excitation 
of double-boundary-localizcd SPPs in the symmetrically- 
sandwiched film. We study in detail the fine structure 
of the two-humped resonance maxima (minima) of the 
transmittance (reflectance), stressing the role of long- 
range and short-range SPPs. As far as we know, the 
fine structure of two-humped resonances corresponding 
to the symmetrically surrounded film has not been stud- 
ied thoroughly in experiments. In the Appendix below 
we draw an analogy between the ELT and the forced os- 
cillations of a well-known classical system of two weakly 
coupled linear damping mechanical oscillators. 



The problem under consideration is of profound in- 
terest not only from the view point of a pure physics, 
but can find numerous applications in optical subwave- 
length devices design 4,50 . For instance, recently, SPP 
in sandwiched structures found itself as surface plasmon 
resonance optical sensors, 51-53 and some of them are al- 
ready commercially available. These sensors are based on 
the interaction between waveguide and SPP modes. The 
basic feature of the sensors is measuring the refractive 
index variation in biological and chemical real-time high- 
definition sensing. It is important to mention that ELT 
phenomena are studied in the microwave (THz) region of 
electromagnetic waves 49,54,55 which attracts a great deal 
of attention. 
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II. PROBLEM STATEMENT AND MAIN 
EQUATIONS 

A. Analytical approach 

Consider an arbitrary polarized plane monochromatic 
wave with wave vector k incident onto a surface of a 
double periodically modulated metal film surrounded by 
dielectric media with permittivities e T , r = ±, from the 
medium corresponding to r = — . Imply that the pe- 
riodicity is caused by modulation of the dielectric per- 
mittivity of a conductor, e — s(r t ), r t — (x,y), so 
that e(r t ) — e (r t + mip 1 + ■m2p 2 ), Pi 2 are elemen- 
tary translation vectors, see Fig. 1. In what follows one 




FIG. 1: (Color online)(a) Geometry of the problem. Diffrac- 
tion by the nonsymmetrically-sandwiched 2D modulated 
metal film (periodical modulation is due to the cylindrical 
inclusions), (b) Fourier spectrum of a periodical array with 
pi = p2 — p, a/ p — 1/3 shown in (a), see Eq. (38). 

has to deal with the Fourier expansion of the function 
£(r t ) = \fi/e{r t ), where e is the mean value of the dielec- 
tric permittivity of the metal, e — (e(r)). The expansion 
of this function over \e — e(r t )]/e coincides in the ze- 
roth orde r with the expansion of the surface impedance, 
V V £ ( r t)j over the same parameter. Therefore, let us 
recall, for brevity, the function £(r t ) surface impedance. 
We write its Fourier representation as follows 



M 



i M exp[i(m x gi + m 2 g 2 )r t ] , (1) 



where £0 = 0, M. is a vector index, or a multiindex, 
A4 = (mi, 1x1-2) (integers mi, m 2 indicate a number of el- 
ementary translations along vectors of reciprocal grating, 
gi, g 2 ); the zero multiindex is O = (0,0). The electric 
fields in the dielectric media are written in the form of 
Rayleigh expansion [the time dependence exp (— iut) is 
omitted everywhere]: 

E r (r) = 5 T -Eexp(ikr) + 

E E^, exp [ik Mt r t + irk T \ Mz {z - S T , + d)] , (2) 



M 

for z > d (z < 0) if r = +(— ). Here d is the film 
thickness, E denotes the electric field amplitude of the 
incident wave, r — (x,y,z), tangential, k t , k^t, and 
normal, k z , fc r |A-fzj components of the wave vectors of 



spatial field harmonics, k 



v k~ , k 



M 



e zk T i_\4 Z , are 



fc_cos#, k t = fc_(sin 0, 0, 0), 
k Mt = k t + migi + m 2 g2, 



oj/c, (3) 



where 6 is an angle of incidence, 
3ie(k T \Mz),^m(k T \Mz) > 0. Analogously to 31-33 
we will take into account the modulation in the bound- 
ary conditions only, so that within the conducting film 
we seek the solution in the form 



E(r) = ^2 E M exp(rfcz + ik Mt r t ), 
k = kV^S, < z < d. 



(4) 



Introducing polarization unit vectors: 
T |+ _ e z x k Mt r |_ e 



M x 



~M 



k\it 
e z x k t 



1+ x k 
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where e^ CT correspond to TM (TE) or p (s) polarization 
for a = — (+) in .Mth diffraction order in the dielectric 
media r, and e _ ' CT are polarization basis vectors for the 
incident wave. With respect to the polarization vectors, 
the electric and magnetic fields in the dielectric media 
are 
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where a = —a; and analogously inside the film. Let 
us introduce the polarization transformation coefficients 
(TCs), T$ a ', as 
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r\a 
M 



E T - 



r\aa L 
M 



E° 



(7) 



Then, excluding from the Maxwell equations and the 
boundary conditions the internal fields, we arrive at the 
following infinite linear system for the TCs corresponding 
to the outer fields 
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M' ,r' ,a' 



r~.r t' 1 00" rpT ' I a" a' 
MM' 1 M' 



V, 



M 



(8) 



The matrix D 



MM' 



and the right-hand side vectors 



V. 



M 



are linear with respect to the modulation, 
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rr'\aa' 
MM' 



^MM'Saa'b 



M 



''MM' 



^m" 7 = &M,obaa'V^ " 



- M 



(9) 



The modulation-independent terms are diagonal both 
with respect to the diffraction order and polarization. 
Explicitly, 
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b M 



= t a 



5rr'£ T ^" Pl( M ti& + {a)' ' £,o(Pt>\mV£^) 1+2 ° (cOSh$)" 



V ' - 



6 T ,-£- +/ /Lj£th$ + ot^t eo(/3-|oV^)^ (cosh*)" 



(10) 



Here 



$ = fed, 



the real part of $ is the film thickness in the skin-depths, 
f3 T \M is the dimensionless z component of the A^th field 
harmonic wave vector, which is related with the core- 
sponding tangential component olm as 



Pt\M — 



k 



V \Mz 



l M 



The nondiagonal terms are 



, a M = k Mt /k. (12) 



iTt' I (7 (j' i 

a MM' — T 



J MM 
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(13) 
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mjvllS^MM') are sines (cosines) of the angle ipMM> 
(a.Mi a M'} between vectors olm, a M'i 

a M ■ at M > . 
=costp M M', & 

e z ■ (ol m x ol M ') 
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where r = +(— ) is for metal-substrate (superstrate) SB 
SPP. The symmetric surrounding (for instance, a free- 
(H) standing film) presents a specific case because of the co- 
incidence of the Eq. (16) solutions for different t. 64 Then 
the initial SPPs existing at the boundary of the metal 
and dielectric half-space become coupled due to the finite 
film thickness, and one obtains the two double-boundary- 
localized (DB) SPP modes: long-range (LR) and short- 
range (SR) SPPs 56,57 . For e + = e_ = e, and, respec- 
tively, /3 + im = /3-\m = Pmi °nc finds from Eq. (15) the 
two roots, @m= Pmi 

/3i, = -^tanh($/2), (3 S M = coth($/2), (17) 

that is the single-boundary localized modes coupled into 
double-boundary localized ones. Therefore, the frequen- 
cies of initial SPPs are repulsed, the spectral degener- 
ation vanishes, and one arrives at two different eigen- 
frequencies described by Eq. (17). Here superscript "I" 
("s") corresponds to the LR (SR) SPP. The LR (SR) 
mode possesses a high (low) frequency and is related to 
-|- 5 antisymmetric (symmetric) with respect to the midplane 

z — d/2 surface charge distribution and the spatial dis- 
= sin^vfA-e, cr = — ■ tribution of the tangential to the film faces electric field 

OLM a M' component. As a result, SR SPP possesses higher Ohmic 

losses. 

Since the modulation is presumed to be small, the 
eigenmodes of the modulated film ( "dressed" modes) dif- 
fer slightly from those existing in the unmodulated film. 
The dispersion relation of the dressed SPP modes de- 
fines the resonance conditions. However for identification 
of the resonance type the modulation can be neglected. 
Bearing in mind that Pm\t depends upon the angle of 
incidence, 9, and the wavelength, A, it is possible to con- 
sider the imaginary part of (16) as the "resonance curve" 
in 9 — A plane. For instance, in the case of rectangular 
symmetry, gi_Lg 2 , gi = .92, imaginary part of Eq. (16) 
reads 

(sin^cosV' + wiKi) 2 + (sin#sin?/; + m2K2) 2 = ^, (18) 

where ip is the angle of the incident plane orientation rel- 
ative to gi (say, "tilting" angle), Ki i2 = gi, 2 /fc = \/pi,2, 

and K T = \J £ T + £?-£e> 2 is the dimensionless wavevec- 
tor of SPP. We designate the curve given by Eq. (18) 
as (7711,7712)1-- All points of 9 — A plane may be classi- 
fied as follows. Indeed, the resonance curves possess the 
finite width due to the dissipative and radiative losses 
and, therefore, it is more appropriate to refer to the res- 
(16) onance vicinity, but not to the resonance point. If a 



(14) 

The resonances in the system are due to existence of 
the eigenmodes in the film, i.e., SPPs. When an evanes- 
cent field harmonic is close to the grazing one, its ampli- 
tude increases significantly as the process of eigenmode 
excitation occurs. For infinitesimal modulation the eigen- 
modes of the film are initial SPPs, with the dispersion 
relation corresponding to the determinant of the matrix 



vanishing in zeroth order approximation in 
modulation. Then det ^^mmv I becomes an infinite 



tt'\o 1 
M 



determinants, so that each of them 



corresponds to some SPP eigenmode of the unmodulated 
film. The equation det CT | =0 has the physical roots 
for a = — only, that corresponds to p polarization: 

((3 + \M tanh $ + £ c ) ((3-\ M tanh $ + £0) cosh 2 $ 

-&=0. (15) 

For a rather thick film, exp($') ^> 1, SPPs existing in the 
film are close to those existing at the boundary between 
the metal and each of the dielectric half-spaces, and are 
single-boundary-localized (SB) SPPs. These modes are 
governed by the dispersion relation 

Pt\m + £0 = 0, 
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point does not belong to any resonance curve, we have 
a nonresonance diffraction case. If a point belongs to a 
single curve with a fixed M. and r, we obtain a single 
diffraction- order single-boundary resonance (SSB). If a 
point corresponds to intersection of several curves having 
different multiindexes {A4, A4', A4", ...}, but the same 
r value, we obtain a multiple diffraction- order single- 
boundary resonance. Note that for specific geometry of 
high symmetry, -0 = and tp = 7r/2, the curves for 
different signs of m 2 and m\ coincide; when ip = ir/4 
curves (m 1 ,m 2 ) T and (m 2 , m i)T are indistinguishable as 
well. Intersection of two resonance curves with differ- 
ent r yields the point of a double-boundary resonance. 
When DB SPPs have a unique multiindex A4 we then 
arrive at a single diffraction- order double-boundary res- 
onance (SDB). SDB resonance occurs only for the sym- 
metric surrounding of the film, e + = e_. Here the ini- 
tial surface modes are coupled mainly through the fi- 
nite film thickness and, in general, the corresponding de- 
pendencies of the rcflcctance/transmittance are of two- 
valley (two-peak) shape due to splitting of LR and SR 
modes. The explicit form of the resonance curve corre- 
sponding to LR (SR) SPP resonances coincides with that 
given by Eq. (18), if one changes £e> to ^tanh($'/2) 
coth($'/ 2 )] in the designation of K T in Eq. (18). 
When DB SPPs arc related to different multiindexes 
(multiple diffraction- order double-boundary resonance), 
which occurs under very specific conditions, they are cou- 
pled through periodicity and the finite film thickness si- 
multaneously. 

Now take up the solution of the system (8). In a 
small region of 9 — A plane, which is of order of the reso- 
nance width, the solution is strongly dependent upon the 
number of resonance curves passing through this region. 
Hence it is convenient to subdivide the set of the diffrac- 
tion orders into a resonance subset, 3?, which contains the 
multiindexes relative to the above-mentioned curves, and 
a nonresonance subset, 01. Thus, we subdivide the initial 
infinite system into the resonance subsystem containing 
the resonance TCs and the nonresonance subsystem con- 
taining the nonresonance TCs. The resonance TCs cor- 
respond to p-components of the TCs with the resonance 
multiindexes, and the nonresonance TCs correspond to 
s-components of amplitudes with the resonance multiin- 
dexes, Tp +<J , and both p- and s-components of TCs with 
nonresonance multiindexes, Tp aa ' , Af em. 

In the main approximation 3 1-33 , which assumes retain- 
ing the quadratic in modulation amplitude terms in the 
matrix elements, and the linear terms in the right-hand 
sides, the resonance subsystem becomes: 



where 

d Tt' r l7 T '\~ I J TT 'I 

u nn' - °nn'O n + a nn , 

jrT"\-<j" /T--l\ T " r '" l<7 " ,r'"r'\a"- fon , 
2^ d KM [ b ) d MW > (2°) 

M,cr" t",t"' 

t',t" 

Here b = the sum with the ovcrlinc, means 

that the terms with the superscript a" = — and a reso- 
nance diffraction order M. = 1Z <G 5i have to be omitted. 
The function S-sro is equal to 1, if within the resonance 
indexes there is the zeroth one, and to otherwise. The 
zero-value of matrix B determinant yields a dispersion 
relation of the SPP modes in the film in the main approx- 
imation. As was discussed above, the block, ' , being 
diagonal relative to the diffraction order, contributes to 
the unperturbed dispersion relation. The nondiagonal, 

linear-in-modulation block, df^P , contains the "inter- 
resonance" or "coupling" modulation harmonic, tji-w, 
which is chiefly responsible for the SPP dispersion curve 
splitting, and for the appearance of the spectral band 
gap 58 ' 59 . The third term in (20) is quadratic in modula- 
tion, and describes the second-order scattering processes 
which provide the broadening and the shift of the disper- 
sion branches. 

The nonresonance TCs, Tp a(7 , are expressed in terms 
of the resonance TCs as 

t',t" 

(22) 

Tp a is a transmission (for r = +) or a reflection (for r = 
— ) coefficient corresponding to p (for a = — ) or s (for a = 
+) polarization in the case of an unmodulated film, s- 
components of the TCs with the resonance multiindexes, 
Tp +a , arc expressed in terms of p-components, Tp a , 
and are small as compared with the latter; they are of 
no interest, and we do not present them herein. The 
diffraction efficiencies are 

|E+(| 2 $te(k +lMz ) |E^| 2 %e(k_\ Mx ) 

™ = ^ kz ' PM = ^f 

From the energy conservation it evidently follows that 

l-J>*f+7Vu]=7>>0, ( 23 ) 
M 

where V ~ / e"|E| 2 dz is the absorbed part of the energy 
flux density. 



B. Single diffraction-order resonance 



^rr'^k'' a ~ ^r <T > (19) -^ a sertain point of A — 9 plane belongs to one or two 

T t^Z! curves (18) having a single multiindex 1Z = (ri,T2) (for 
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other parameters being fixed) the vicinity of this point 
defines a single diffraction-order resonance. Moreover, 
according to the above classification, in the case of one 
curve (for a single t), we have a SB resonance, and if 
the point is the intersection of two curves with different 
r values, we have a DB resonance. The resonance TCs 
are then similar to those obtained for ID modulation, cf. 
Ref. 33, 



" T t\-+~ 
1Z 




rpT\ 





cos 9 sin ipizo 
cos Ipno 



L^inieM-n^ , (24) 



where 



^1Z 



L n = 2r(^, K - Sr.+T^) /A n , 
Pr\n = Pr\n tanh $ + & + G T n , 
/3 + | TC /3_| TC -T^cosh- 2 $, T 



iz 



g tz 



(25) 
(26) 



G iz > 
(27) 



GT 



K 



E 

AT 



£,k-n£,N-iz ( 2 , 

COS IpTZM- 



Er/^l^sin 2 Ip-RAT)- 



Stress that T n l+a 



(28) 

0{£k), \Tt°\ » The 
resonance TCs have two poles. For a nonsymmetric sur- 
rounding, these poles are releted to SPPs existing at the 
opposite film faces, while in the case of a symmetric sur- 
rounding, they are related to LR and SR SPP modes. 

The structure of the resonance TCs indicates that the 
coupling strength between the incident wave and the SPP 
excited is proportional to the scalar product of the SPP 
and the incident wave magnetic fields. Note that the 
SPP magnetic field is orthogonal to its propagation di- 
rection and is parallel to the interface so that the SPP 
is effectively excited by the projection of the tangential 
component of incident wave magnetic field onto H^ t (or, 
alternatively, by the scalar product of tangential compo- 
nents of the electric SPP fields and the incident wave) . As 
for p(s)-polarization of the incident light, the projection 
is H t cos ipfco (Ht cos 9 sin tpno), where ipno is the angle 
between K nt and H t . For instance, with the incidence 
of purely p- (s-) polarized light, the SPPs propagating 
parallel (perpendicular) to the incidence plane cannot be 
excited, even if wavevector k-jit is close to the poles in 
Eq. (24). 

The zeroth-order polarization matrix is 



'rpT\ + + rpT\ + - 
rpT\ h rpT\ — 



-F 








Tl 



2 cos 9 T cos 9 sin 2 tpno — cos 9 T sin 2ip-jzo 
—Tcos9sin2ip-ji0 2t cos 2 ipno 



^0|TC'( 29 ) 



where 



j O\TZ 



IZ 



cos 9 t Atz 



Pt\iz + (Pt\K 



T K J (cosh$) 7 
x (cosh$) _J ^ 



• (30) 



9 T is the angle of propagation of the zeroth-order wave 
in rth dielectric media relative to Oz axis (in the super- 
state 9- = 9). The structure of these coefficients shows 
an interference caused by the competition between the 
nonresonance channel (the terms Tp a ), and the reso- 
nance channel (the second terms in Eq. (29), which is 
proportional to ^iz^-tz/^tz)- The nonresonance term 
should be retained for the zeroth-order reflectance, since 
it is of order of unity, |T F ' CT | ~ 1, while for a rather 
strong resonance it can be neglected for the transmit- 
tance in the vicinity of the resonance maxima, since 

I^f' CT | ~ l£o| ex P( — *C 1 and is much smaller than 
the resonance input into the transmittance. 

The wavelength resonance width, AA/A, is contributed 
both from the dissipation losses, being proportional to 
and from the radiation losses due to SPP scatter- 
ing into the outgoing propagating waves. It is of or- 
der AA/A ~ Wo\Wo + C(|| 2 |)], as it follows from (24)- 
(30). The term proportional to 0(|£ 2 |) is caused by 
the radiation losses and in the simplest case may be 
represented approximately as J2m QaT&^-AT&V'-tc/A-IAA 
where |CV| ~ 1, and the summation is performed over 
those Af that satisfy $sm((3 T \ss) = 0. According to 
the formulas describing the resonance and zeroth-order 
TCs, the optimal amplitude of the resonance harmonic 
is \£n\ ~ v^o an d this leads to the resonance width 
of order AA/A ~ Col^ol- The optimal amplitude is re- 
lated to the maximal SPP excitation and, consequently, 
to minimal/maximal value of reflectance/transmittance. 
In spite of the fact that we assume the modulation to be 
rather small, let us make a rough astimation of the reso- 
nance width for the hole arrays. For Ag films in the visi- 
ble and near infra-red frequency region the impedance of 
the film is ~ 10 _1 , while for the holes the impedance 
|&| ~ 1. Therefore, \£ n \ ~ |A£| ~ 1 » and 
AA/A ~ |^| ~ 10 _1 . For wavelength of order A ~ lfj,m 
the resonance width can be estimated as AA ~ lOOnm, 
which is in good qualitative agreement with numerous 
experimental results. In the experiments the resonance 
width is equally affected by the non-plane character of an 
input light wave and the finitencss of the periodic array. 

Note that the resonance width is very important for 
problems of sub-diffraction-limitcd optical imaging (op- 
tics of volume or surface "superlenses"). Sufficiently 
broadened SPP resonances may be efficiently used to 
enhance the evanescent modes and thus to recover the 
subwavelength information on nanoobjects, see Ref. 61. 



C. Multiple resonances 

The approach developed allows one to consider the 
diffraction problem for the resonances of arbitrary mul- 
tiplicity on 2D-periodical structures possessing an arbi- 
trary symmetry. But from the experimental point of view 
the resonances of fourfold multiplicity, which take place 
for the normal incidence, 9 — 0, at the square and rectan- 
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gular periodical arrays are of special interest. Here, and 
in what follows, we will concentrate on the structures 
possessing C 2v symmetry: £(C 2v r t ) = £( r t)> the geomet- 
rical symmetry (Brillouin zone symmetry) is supposed to 
be Ci v , i.e. gi is perpendicular to g 2 and their modules 
are equal, gi = g2. As it seen from Eq. (18), at normal 
incidence the multiple S PP reso nances arise at the wave- 
lengths A£ x r2 = pK T j ' \Jr\ + rf, where p is the period of 
the structure. Points 9 = 0, A = A£ r . 2 in the — A plane 
are the intersections of four resonance curves for the sin- 
gle r (SB resonance) or eight curves for both r = ± (DB 
resonance in symmetric surrounding). We mark a four- 
fold SB resonance as [ri,r 2 ] r , which corresponds to the 
intersection of the following resonance curves: (ri,r 2 ) T , 
{-ri,r 2 ) T , {n,-r 2 ) T , and (-ri, -r 2 ) r . A fourfold DB 
resonance is marked as [ri,r 2 ] respectively. In addition, 
if all values of a certain function F_m with subscripts from 
the above subsets are equal, we denote it by i*j rijT . 2 i. 

Similarly to the single resonance, when we have the 
fourfold one, the pair of field harmonics with wavevectors 
k (rur 2 )t = rigi + r 2 g 2 and k ( _ n _ r2 ) t = -ngi - r 2 g 2 
is efficiently generated by the projection of the incident 
wave magnetic field onto the direction perpendicular to 
k( ri)T . 2 ) t (the direction, corresponding to the SPP mag- 
netic field with this wavevector), viz. by H sm(j) ri ^ 27 
where <f>-jz is the angle between k-jn and H 4 = H. Analo- 
gously, the amplitudes of the field harmonics correspond- 
ing to diffraction indexes (±n,=Fr 2 ) are proportional to 
H sin 4>- ri t r 2 . Thus, for special polarization of the in- 
cident wave, at cpn — 0, ±7r, ±lZth the resonance field 
harmonics are not excited via the first-order scattering 
process (but rather via the higher-order processes). This 
reduces the multiplicity of the resonance twice in the 
main approximation. For an arbitrary polarized inci- 
dent wave, the polarizations of the zeroth-order trans- 
mitted and reflected waves are formed mainly by the in- 
terference contribution of all resonance field harmonics 
due to single back-scattering. Also, they are contributed 
from the zeroth-order components related to the "scatter- 
ingless" reflection and transmission for an unmodulated 
film. Otherwise, the polarizations of both propagating 
and evanescent nonresonance field harmonics are formed 
mainly by the single scattering from the zeroth diffrac- 
tion order and from all the resonance diffraction orders. 

For the structures possessing C 2v modulation symme- 
try, the solution of (19) is considerably simplified in the 
vicinity of normal incidence. In the case of [r, 0] or 
[r, 0] T resonance, the resonance TCs are similar to those 
given in Eq. (24), cf. explicit expressions in Appendix B, 
Eq. (Bl). Though in f3 T \n a linear-in-modulation term, 

£ 2 K, arises, which is the inter-resonance harmonic respon- 
sible for the splitting and shifting of the resonance. 

To get a better understanding of the resonance diffrac- 
tion, consider the eigenmodes of C 2v structures. They 
are defined approximately by relation BT = 0. In par- 
ticular, the cigcnfrcqucncics may be found from equation 
det B = 0. For SB-localized SPPs when we have a rather 



thick film, exp(— $') <C 1, the structure of eigenmodes 
can be obtained in the approximation of the half-space 
problem, see the detailed analytical treatment in Ref. 58. 
Thus, we restrict ourselves to the eigenmodes that are 
close to those corresponding to the metal-dielectric in- 
terface and not coupled through film thickness. Bear- 
ing in mind the homogeneous problem setup, we must 
change the notations so that one of the resonance k- 
vectors in the diffraction problem, k^, become the SPP 
quasi- wavevector, q, ending in some Brillouin zone. The 
other wavevectors close to the resonance conditions corre- 
spond to the "resonance satellites" , that compose a cou- 
pled SPP state. 

Concentrate first on the simplest case of a twofold cou- 
pling through the periodicity. Suppose that the reso- 
nance wavevectors are k# t — ► q and k#/ t — > q'. Then we 
can consider the diagonal in r homogeneous subsystems 
of the resonance system (19) which contains the two TCs, 
for the diffraction orders 1Z = (ri,r 2 ) and TV = (r[,r 2 ). 
Using Eqs. (10)-(14) (and assuming tanh <f> = 1), we write 
it as 

( - ^ + f &f C0S ^') ($>)= 0. (31) 

Here /S.1Z = 1Z — 1Z' , and we neglect the quadratic- in- 
modulation amplitude terms, supposing that they do not 
exceed the linear-in-coupling harmonic £,±ak term. We 
take into account that in the rth dielectric half-space the 
SPP electric field has predominantly a z-component and 
therefore, change a to E^ z . From this system we 
see immediately that the initial SPPs with wave vectors 
oriented at ip-jzw — ±tt/2 are not coupled, and it is thus 
the case where it is necessary to retain the quadratic 
coupling terms 65 . Note that, as we see, here the 2D 
problem is reduced to the ID one, corresponding to the 
twofold SPP coupling through the harmonic grating of 
the period 2n/\kn lt — k^ 2t | (and with amplitude £±ak), 
where the quasi-wavevector of the dressed SPP, q, and 
the satellite, q' = q+ (r[ — ri)g± + (r' 2 — r 2 )g 2 , refer to the 
sides of a Brillouin zone, see Fig. 2. Then fl T \ n = [i T \ n i = 

(3 T , where (3 T — e T — (qc/uj) 2 /e T , and the dispersion 
relation becomes 

[3 T {UJ ± ) = -fa T ^|Arl-Ar • | cos tpnw I ■ (32) 

Remind that ip-R.il' is the angle between q and q'. Hence, 
we obtain the two roots: the first one is for the high 
frequency mode, u + , and the second one is for the low- 
frequency mode, uj~ . Specifically, neglecting the modu- 
lation of the small real part of the surface impedance 
as compared with that of its imaginary part we have 

|-m = -£m> and \J i A ni-An = «||atc|- Hence, the 
eigenfrequencies are 

w ± /w p h ^ 1 + e T S,a/2±ie T S,o\^ATzcosipTzn'\, 

u ph = qc/y/e^ , (33) 
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FIG. 2: (Color online) An example of a twofold coupling of 
the SPPs in the case where A1Z = (1, 1). The upper figure is 
for the 27? problem, and the lower is for the similar ID one. 

where the upper (lower) sign is for high (low) frequency 
SPP. It is important to note that the quality(Q)-factor 
for the high-frequency mode is higher than that for the 
low-frequency mode. 

For the specific choice of coordinate origin, such that 
£ak = *|£aw|, f (r t ) = Co + + 2i\£ An \ cos[(q - q')r] + 
the electric field amplitudes of the eigenmodes obey 
the relation (ET, n / E^,)^ = ±sign[cosr/>7i7i']- As the 
coupled SPPs propagate in the opposite directions, ip = 
±7r, we have q = q', and the field structure at the 
interface for the eigenmodes is 

£ z T (-)(r t )~cos(qr), E<+\v t ) ~ sin(qr), (34) 

Note that under the condition of the prevalence of the 
coupling in the first order scattering, or, in other words, 
when the coupling harmonic dominates over the others, 
\S,atz\ > \£m\ 2 : the structure of the spatial field distribu- 
tion is given by the coupling harmonic. In other words, 
the field maxima of the high frequency mode coincide 
with the "less metallic" regions (where |£| is higher) rel- 
ative to the coupling harmonic and vice versa for the low 
frequency mode, see Fig. 3. 

Discuss the eigenmodes for fourfold SB SPP resonance. 
We limit ourselves to [r, 0] resonance, such that k( r o)t = 

-k(- r ,0)t; k(0,r)t = _ k(0,-r)t, and k( 0i±r ) t _L k(± r) o)f. 

Remind that the resonance amplitudes corresponding to 
perpendicularly oriented vectors are not coupled in the 
first scattering order. One can make sure that four eigen- 



FIG. 3: (Color online) Spatial distribution of the squared elec- 
tric field z-components of the eigenmodes. The solid (dashed) 
curve is for the low (high)- frequency mode. The "more metal- 
lic" regions of the coupling harmonic contribution are shown 
by the darker regions of the contourplot. Pluses and minuses 
refer to the signs of the surface charge density. 

modes exist. Two of them are "mixed", that is, all res- 
onance amplitudes are nonzero: EJ, r Q . , ET ±r s z ^ 0, 
moreover, ET „s = ET „s , and ET n s = ET n s . 

' (r,0)z (— r,0)z> (0,^)2 (0,— r)z 

Thus, the field structure has the form of the linear com- 
bination of two cosines having their periods along gi and 
g2- These modes cannot be excited at normal incidence. 
Other two modes, which may be excited, have zero res- 
onance amplitudes: for one mode ET rQ ^ z = — ET_ r ^ , 
EJ n ,, = 0, and ET n , = -ET n , ' ET, m = 6 for 

(0,±r)z ' (0,r)z (0, — r)zi (±r,0)z 

another one. Consequently, the spatial structure corre- 
sponds to sinuses. Thus, the structure of the eigenmodes 
in this particular instance of fourfold SB SPP resonance 
is similar to that of twofold SB SPP resonances; they 
have the form of a pair of standing SPPs, that corre- 
spond to the sinus-type spatial field distribution, and 
higher- frequency (and also higher Q-factor) branches. 
This means that in Eqs. (33) taking "+" we should set 
| cos tprni' I - ¥ lj Ar — > (2r, 0), and the first (second) 
standing SPP quasi- wavevector, q, becomes parallel to 

k(r,0)t [k(o,r)t]- 

The solution of the inhomogeneous problem is the 
superposition of the above standing SPPs with the 
"weights" proportional to cosines of the angles between 
HT. w, HT Q w, and H. The eigenmodes for other SB- 
fourfold resonances may be treated in a similar way, they 
are the combinations of simplest SPP modes as well. 

For DB-localized dressed SPPs the field structure is 
defined by coupling of the initial SPPs both through 
the modulation and through the finite film thickness. 
The field structure resulting from the coupling through 
the film thickness may be understood from the exam- 
ple of the undressed SPP existing in the unmodulated 
film. The eigenfrequencies for this case are defined by 
Eq. (17). Note that inside the metal the amplitude of the 
tangential-to-interface component of the electric field is 
higher than z-component, \Et/E z \ ~ \y/e\ 3> 1, in con- 
trast to the dielectric half-spaces. Remind also that in- 
dependence of the electric field inside the film for LR and 
SR eigenmodes is 1,2 

E { t L) ~ sinh(fcz - $/2), E { t S) ~ cosh(kz - $/2), (35) 
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where (L) stands for LR SPP and (S) for SR SPP. 

For dressed DB-localized SPPs with the twofold cou- 
pling, the dispersion relation is 



P (L±) = -(fo±£ AW |cos^|)tanh($/2), 
/3 (S±) = -(& ± |ak| cos ^|) coth(*/2), 



(36) 



where ft-) = - (qc/uj) 2 /e d , e+ = e_ = e d , su- 
perscripts in brackets mark the eigenmodes. Here the 
initial SPPs (and corresponding diffraction indexes) are 
the same as in the above case of the dressed twofold SB 
SPP. The spatial field distribution for eigenmodes inside 
the film for ip = it (q = — q') has the form 



1T 



Sf->(r) 



(r) ~ e q cosh(fcz — $/2) sin(qr) 
e q sinh(fcz — $/2) sin(qr) 



(r) ~ e q cosh(fcz — $/2) cos(qr) 



E| i+) (r) ~ e q sinh(£;z - $/2)cos(qr) 



(37) 



z-components of the electric field in the dielectrics at 
the film faces are of the same form as in Eq. (34), i.e. 
££ (s-) (r) ~ El {L ^{r) ~ cos(qr) and E T z {S+) (r) ~ 
El (r) ~ sin(qr). It is possible to obtain the eigenfre- 
quencies and the field structure for the fourfold DB SPP, 
proceeding in the same way, in which we generalized the 
twofold SB SPP to the fourfold SB SPP. 

Analyzing the denominator of TCs, on can make sure 
that the excited DB SPPs are related to "+" modes only. 
Indeed, with the film thickness tending to infinity, we can 
see that the two poles of TCs in Eq. (Bl) of Appendix B, 
A n = 0, become /3 T \ n + £e> + £,2iz = 0, in neglecting 
the second-order scattering processes. Choosing an ori- 
gin such that £27?, = ^I&tcIj the poles coincides with "+" 
branch of Eq. (32), where we have instead of £ah, 
and I cos ipnw 1 = 1- This corresponds to the upper fre- 
quency branch, Eq. (33), and, hence, to the "sine" field 
distribution in the x — y plane. 

Under the assumption of C^ v modulation symmetry 
the zeroth-order TCs for SB [r, 0] T resonance have simi- 
lar to Eq. (29) structure as well (see Appendix B). How- 
ever, the difference consists in identity T^ aa — 0, which 
ensures that the polarization of the zeroth-order trans- 
mitted wave coincides with that of the incident wave. 

Along with the anomalies related to SPP excitation 
there exist Rayleigh anomalies. They correspond to the 
boundary between homogeneous (propagating) and in- 
homogeneous (evanescent) waves in different diffraction 
orders, viz to the vanishing of z-component of one of the 
wave-vectors, k T \ Mz = 0, or 0j^i T — 0. From the math- 
ematical point of view, the Rayleigh anomalies are the 
branch points; they set the conditions for the breaking of 
the incident angle or wavelength derivative of the trans- 
formation coefficients. In what follows we will mark the 
Rayleigh anomalies as (m^TT^)^ (which corresponds to 
the vanishing of (3 T \ miim2 ) or as [m 1 ,m 2 ]^ ? (which corre- 
sponds to the vanishing of l3 r \[ mi ,m 2 ])- 
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FIG. 4: (Color online) Modulation harmonics £ mi , m2 depen- 
dence upon a ratio of the cylindrical inclusion radius to the 
period of the square array, a/ p. The arrows 1, 2, 3, 4 indi- 
cate a/p values of the hole arrays studied in Refs. 5, 6, 14, 9 
correspondingly. 



Since modulation harmonics play a crucial role in SPP 
excitation, et us consider below the Fourier representa- 
tion of the structures predominantly used in the experi- 
ments. 



D. The modulation spectra 

The majority of experimental works deal with hole ar- 
rays fabricated with the use of metal films deposited onto 
a dielectric substrate (predominantly onto quartz), see 
Fig. 1. Instead of holes we will consider a cylindrical 
inclusions in the film. The inclusions are supposed to 
consist of a metal or a semiconductor with the dielectric 
permittivity different from the permittivity of the film. 
Both the inclusions and the film must be highly conduct- 
ing. This structure, as we will see below, may qualita- 
tively describe the optical properties of nanohole arrays, 
even thought the film does not contain holes as such. 
Furthermore, the results from our studies may be used 
for modeling and describing arrays of cylindric nanopar- 
ticles, so far examined experimentally. Indeed, the struc- 
ture considered may be fabricated by making an array 
of cylindric nanoparticles with impedance which are 
immerged into a conducting film. 

For the inclusions of a round cross-section with radius 
a, we have 

- 2^A£ ira 2 At; 
£m = ^ Ji(agM), 4o=C/H ^ — , 

Sm = migi + m 2 g 2 , £ = \p t x p 2 \, (38) 

where A£ = — is the difference between film (£/), 
and inclusion (£,) surface impedance, J\ is the first-order 
Bessel function. A typical example for the Fourier spec- 
trum of the array of cylindrical inclusions for the square 
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symmetry structure is shown in Fig 4. The square sym- 
metry structure means, that the angle between transla- 
tion vectors is 7r/2 and their modules are equal, p\ = P2- 
In the following calculations we take to be equal to the 
impedance of Ag (using the wavelength dependence from 
Ref. 62); the impedance of inclusions will be modeled as 
& = w £/i where w is the dimcnsionless parameter. 

As evident from Eq. (24), the amplitude of the res- 
onance (IZth) diffraction order is proportional to 
H^ a oc therefore, the efficiency of excitation of SPP 
in this order is strongly dependent upon the correspond- 
ing "resonance" amplitude For large values of ri,r% 
the Fourier amplitude tends to zero as 



-3/2 



\n\ > i. 



(39) 



Thus, the resonances in high diffraction orders are less 
efficient. If the inclusion diameter is far smaller than 
the modulation periods, the low-order amplitudes are ap- 
proximately independent of their order: 

£m ~ 7rACo 2 E- 1 [l + 0{ag M )} for ag M < 1. (40) 

So, for very thin inclusions the Fourier amplitudes of the 
structure show a slight decrease, which is valid for an 
arbitrary cross-section of the inclusions. In its turn, it 
significantly broadens the SPP resonances: they become 
"diffusive" due to different scattering processes of the 
excited SPPs. 

Let us write also the Fourier expansion of the structure 
having the inclusions of a rectangular cross-section, 



4aia 2 A£ 



(»i • gx)(&2 • gjw)S 



. ai • g M . a 2 • g M 
sin sin . 



(41) 



where ai,2 are the vectors oriented along the sides of 
the rectangular inclusions with modulus equal to their 
lengthes, ai 2- 



III. SB SPP RESONANCES. NONSYMMETRIC 
SURROUNDING 

Bearing in mind the experiments already made, 
periodically-modulated films surrounded by two different 
dielectric media are of special interest. These structures 
can support both SB-localized SPPs and, under the spe- 
cial conditions, DB-localized SPPs. The latter will be 
discussed separately. In the vicinity of normal incidence 
we mainly focus on, the resonances basically correspond 
to SB-localized SPPs. An example of the transmittance 
and reflectance wavelength spectra for the strictly normal 
incidence is illustrated in Fig. 5. The calculations are for 
Cn v array of inclusions with a round cross-section in the 
Ag film surrounded by air superstrate and quartz sub- 
strate. The transmittance and reflectance extremes are 
pronounced at wavelengths for the SPPs existing on un- 
modulated metal-quarts and metal-air interfaces. Note 
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FIG. 5: (Color online) Dependence of transmittance and re- 
flectance upon the wavelength for the wave normally-incident 
from the air superstrate onto the silver film bounding with 
the quartz substrate. The parameters of Ca v array, are 
pi = p2 = p — 0.9pm, a/p — 1/3. The incident wave electric 
field E is parallel to gi (to Ox axis, see. Fig. 1). The inclusion 
impedance, £j, is taken to be £i = 2£/. The film thickness is 
equal to 2.3 skin depths. The transmittance and reflectance 
of the unmodulated film are plotted with dotted lines. The 
vertical solid and dashed lines indicate the wavelengths rele- 
vant to SPPs on the unmodulated boundaries and Rayleigh 
points respectively. 



that while the extremes are shifted with respect to ini- 
tial SPPs wavelength, wich is due to the modulation and 
finite film thickness influence, Rayleigh anomalies are un- 
shifted. 

Now discuss the vicinity of the strictly normal inci- 
dence. In Figs. 6 (a-c) and (e) the transmittance contour 
plot is shown as a function of the wavelength or/and 
photon energy and the incident angle for different polar- 
ization of the incident wave. In Fig. 6 (d,f) the resonance 
curves (18) corresponding to the contour plots in Fig. 6 
(a-c) are presented. There is an excellent agreement 
between the transmittance features and the resonance 
curves. It should be stressed that the short-wavelength 
region corresponds to strong dispersion, £(A). In Figs. 6 
(a-c) the special symmetry instances are shown. Viz., 
orientation of the incident plane relative to one of the 
reciprocal grating vectors (in our case it is gi) at an- 
gles (or 7r/2) and ir/4 results in the coincidence of 
the resonance curves, see discussion below formula (18). 
The coincidence is appreciable in the vicinity of [1,0] + , 
[1,1]+ and [1,0] _ resonances. In Fig. 6 these vicinities 
are marked by squares. We see that in Figs. 6 (a-c) from 
one to three resonance features belong to these regions. 
Evidently, in Figs. 6 (d,f) from one to three resonance 
curves resonance intersect these regions. Alternatively, 
four resonance "mountain ridges" intersect in the vicin- 
ity of a single point, see Fig. 6 (e), which displays the 
transmittance contour plot for the nonspecific in the an- 
gle tp case. 

Emphasize, that the intersections of solid and dashed 
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FIG. 6: (Color online) Dependence of transmittance upon the wavelength (energy) and the incident angle for different polar- 
ization of the incident wave and different orientations of the incident plane. The array and the film parameters are the same 
as in Fig. 5. In (d, f) solid and dashed lines are the resonance curves, Eq. (18). Solid (dashed) squares indicate the vicinity of 
the quartz-silver (air-silver) resonances which are discussed in detail in the text. 



resonance curves in Figs. 6 (d,f) and the intersections of 
corresponding resonance "mountain ridges" in the trans- 
mittance contour plots are due to DB SPP resonances. 
The latter comply with the excitation of dressed SPPs, 
resulting from the coupling of initial SPPs of different 
film faces through the modulation. 



In the following subsections we examine in detail some 
of SB fourfold resonances. 



[1,0] 



resonance 



Consider first the longest-wavelength resonance (at 1.4 
fim in Fig. 5) in the vicinity of close-to-normal incidence. 
This resonance arises due to excitation of the fourfold 
[1,0]+ SB SPP at the metal-quartz interface having G± v 
modulation symmetry. 

In Fig. 7 (a, c) the enlarged fragments of Figs. 5 and 
6 (a) are shown in the vicinity of [1,0]+ resonance for 
the incident plane orientated parallel to gi and the po- 
larization of the incident wave such that Ellgi. Note 
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FIG. 7: (Color online) (a) The dependence of the zeroth-order transmittance, To, upon the wavelength and incident angle in 
the vicinity of [1,0]+ resonance on silver array of inclusions with the symmetry Cav for E||gi. In (b) the reciprocal grating 
is shown, where the radii of solid and dashed circles are equal to the SPP wave vector modulus in quartz and air half-spaces, 
respectively; the circles separate the regions of propagating (filled area) and evanescent waves in the superstrate and substrate. 
In (c) the wavelength dependence for the zeroth-order transmittance, to, and reflectance, po, are shown, which correspond to 
the grey section in (a). The inset shows the vicinity of the Rayleigh anomaly. The parameters of the array are the same as in 
Fig. 5. 



that the resonance features are heavily dependent upon 
the polarization of the incident wave and the orienta- 
tion of the incident plane, cf. transmittance dependen- 
cies shown in Fig. 6 (a), (b), (c) and (e). For instance, 
for p polarization, when E belongs to the incident plane 
parallel to gi (ip = 0), in the vicinity of normal inci- 
dence, only the pair of the resonance waves are excited, 
see Subsection II C. These waves correspond to (±1,0) 

diffraction orders (TCs are equal to zero). In 

Fig. 6 (a) and in Fig. 7 (a) one can observe the inter- 
section of the two resonance features, relevant to the in- 
tersection of two resonance curves: (1, 0)+ and (—1, 0)+, 
Fig. 6 (d). In deflecting from the normal incidence, the 
projection of H t onto H^ ±1 becomes nonzero, and it 
appears as an feature in Fig. 6 (a) close to the coin- 
ciding resonance curves (0, ±1)+ in Fig. 6 (d). When 
E is perpendicular to the incident plane parallel to gi 
(ip = 0), [s-polarization, see Fig. 6 (b)], the only ex- 
cited pair of resonance waves comply with the diffraction 
orders (0, ±1) (resonance coefficients T^ ±1 ). If the inci- 
dent plane is oriented at ip ^ 0, 7r/4, then two pairs of res- 
onance waves are excited both for p- and s-polarizations 
of the incident wave in the vicinity of normal incidence 
[in Fig. 6 (c,e) s-polarization case is shown]. The reso- 
nance waves correspond to diffraction orders (±1,0) and 
(0, ±1). In Fig. 6 (f) this appears as the intersection of 
four resonance curves: (±1, 0) + and (0, ±1)+, coinciding 
for ip = 7r/4. 

The wavelength dependence of the zeroth-order wave 
has the typical Fano-type profile which consists of the 
neighbouring minima and maxima, see Fig. 7 (c). This 
is due to the interference (see Subsection II B) be- 
tween nonresonance and resonance transmittance mech- 
anisms. The maxima of the Fano profile is red-shifted 
as compared with the wavelength of the SPP existing at 
the nonmodulated metal-quartz interface. This shift is 



mainly due to the scattering by modulation through the 
diffracted inhomogeneous waves, and partially due to the 
finite film thickness. 

The amplitude of the excited resonance waves is pro- 
portional to £[ 10 ] Fourier harmonic amplitude, while the 
zeroth-order transmittance and reflectance are propor- 
tional to its squared value, as it follows from Eqs. (24)- 
(30) (see also Appendix B). Therefore, this harmonic 
plays the most important role in the excitation of the 
SPP eigenmode for [1,0]+ resonance. Actually, the struc- 
ture of the eigenmodes in the modulated film (dressed 
SPPs) define the distribution of the near-field completely. 
As indicated in Subsection II C, the only excited eigen- 
modes are high-frequency ones. For instance, when E is 
oriented along gi, the distribution of the squared field 
on the quartz side take the form |E + | 2 ~ sin 2 [k^ ) 4 r] = 
sin 2 (gir), similary to the spatial field distribution for 
the twofold coupling given by Eq. (34). Indeed, this is 
in compliance with the left-hand part of Fig. 8 (a). In 
the general case, when E is oriented relative to gi at an 
arbitrary angle ip, the near-field structure is formed by 
the interference of the eigenmodes. One of them corre- 
sponds to sin(gir) and the other is for sin(g2r). As a 
result the field on the interface has the form |E + | 2 ~ 
[Xisin(gir) + % 2 sin(g 2 r)] 2 , where the "weight" coeffi- 
cients xi (X2) are proportional to cost/' {simp) respec- 
tively. This interference is clearly seen in the right-hand 
part of Fig. 8 (a). Thus the near- field structure is in 
excellent agreement with the experiment, cf. Ref. 14. 

Consider the field distribution along z axis, see Fig. 8 
(b). Since the far-side SB SPP resonances provide small 
reflectance dips [the lowest reflectance is of order of 93%, 
see Fig. 7 (c)], we notice a clear interference pattern 
along z axis in the air half-space. The field penetrat- 
ing into the film decays exponentially, and excites a SB 
SPP at the quartz-metal interface. The SPP field has 
a well-pronounced "torch" structure. It decays expo- 
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FIG. 8: (Color online) The squared field distribution is shown 
for [1,0]+ SPP resonance at A = 1.4023 fim, 9 = 0. The 
inclusions are shown by the dashed lines. The parameters of 
the array are the same as in Fig. 5. 

nentially into the film at a distance of order of skin- 
depth, S, where 5/ A ~ l/|Vt| <C 1 and decays into the 
dielectric medium (quartz) at a distance of (5+, where 
S T /X ~ ■\/\e/e T \/2w ~ 1. Note that the field pattern 
along x axis in the film is shifted to a half-period with 
respect to the pattern in the quartz. This is because 
the SPP electric field in the metal has predominantly 
tangential-to- interface component as compared with a di- 
electric where the electric field has predominantly normal 
to the interface component. The intensity of the far field 
in the quartz half-space is constant, since there exists 
only one propagating wave (at zeroth diffraction order). 

B. Polarization properties of [1,0]+ resonance: 
dependence upon the inclusion shape 

In this subsection we give a theoretical interpretation 
of the recent experiments 39,41 concerning the influence 
of the hole shape on ELT through periodic hole arrays 
within the framework of our theory. Moreover, we de- 



scribe the transmittance behavior while changing the po- 
larization of the incident light. 

Consider a silver film sandwiched between an air super- 
state and quartz substrate with C% v modulation symme- 
try caused by the inclusions of rectangular cross-section. 
This symmetry reduces to C^ v symmetry for a circle 
cross-section of inclusions. Now let us give a close con- 
sideration to the [1,0]+ SB fourfold SPP resonance at 
strictly normal incidence. As was discussed above, for 
C2v symmetry structures the polarization of the trans- 
mitted radiation depends upon the polarization of the 
incident wave. Suppose that the direction of gi is parallel 
to the small side of the rectangular, a\. If the incident- 
wave electric field is oriented along gi [g2], the incident 
wave magnetic field is then matched with that of (±1, 0)+ 
[(0,±1)+] SPP pairs having wavevectors ±gi [±g 2 ] for 
normal incidence, see Fig. 7 (b). According to our ap- 
proach, the transmittance coefficient, to, is proportional 
to £± x o (£o±i) with a of polarization along gi (g 2 ). 
This leads to t ~ (A£) 4 (a/ p) 4 J^(2na/ p) for round 
cross-section and to ~ (A£) 4 (a 2 /p) 4 sin (wai/p) [to ~ 
(A£) 4 (ai/p) 4 sin (wa^/p)] for rectangular cross-section 
with E orientation along gi (g 2 ) [see Eqs. (38), (41)]. We 
take the same geometrical parameters for periodic arrays 
as in Ref. 39, i.e. p = 425 nm; a\ x a 2 = 75 x 225 nm 2 for 
small rectangles and a\ x a 2 = 150 x 225 ran 2 for large 
ones; and 2a = 190 nm for the circle diameter which 
corresponds to the resonance position at A ~ 0.72 pm. 
Taking the film thickness to be equal to 3 skin depths, 
= 3, and the inclusion impedance to be & = 2£/, 
where £/ is taken for the silver, we obtain the transmit- 
tance of order 2%, which is nearly 7 times higher than 
the transmittance through the unmodulated film. 

When E is directed along gi, the zeroth-order trans- 
mittance amplitude is proportional to a squared £±i,o 
harmonic amplitude. This amplitude increase for large 
rectangles as compared with the circles and this results 
in redshift and in a rise of the transmittance maximum, 
see Fig. 9 (a). This tendency is in accordance with the 
experiment. But for the small rectangles the resonance 
harmonic amplitude is the smallest one and this leads to 
blueshift and a decrease in the transmittance maxima as 
compared with the circles and large rectangles. This is 
not what was observed experimentally. We may elimi- 
nate this discrepancy only by adjusting the impedance 
of the inclusion. We use this adjustment because, in our 
approximation, we do not allow for the influence of the 
modulation upon the volume of the film, but the modu- 
lation appears in the boundary conditions only. In other 
words, the fields decay into the film with the same decre- 
ment at any cross-section by a plane parallel to axis Oz. 
On the other hand, it is evident that under a strong mod- 
ulation, with the impedances of inclusions and the film 
differing greatly, the fields in the inclusion and in the vol- 
ume of the film decay in a different way. For instance, 
inside the holes the field show a weaker decay than in the 
metal regions. Thus, adjusting the impedance, we have 
found that for £j > 2.67£/ the transmittance maximum 
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FIG. 9: (Color online) The wavelength dependence of trans- 
mittance for the inclusion with different cross-sections at 
strictly normal incidence onto Ag film deposited onto the 
quartz substrate. The polarization of the incoming light, E, 
is perpendicular to the long axis of the rectangles in (a) and 
parallel to the long axis of the rectangles in (b). The dotted 
curves are for the large rectangular inclusions. The dashed 
and dash-dotted curves are for the small rectangular inclu- 
sions. The solid curves correspond to the cylindrical inclu- 
sions. The film thickness is equal to 3 skin depths. In the 
insets the amplitude of the (±1,0) and (0, ±1) modulation 
harmonics are shown. The dashed rectangulars in the insets 
indicate the modulation amplitude for the dashed curves. 



discovered polarization dependence of the transmittance 
upon the the ratio of the rectangular holes to a result 
from interaction between SPP resonances at the surface 
of the metal and shape resonances inside the holes. 

We hold to a different, simpler point of view. It is ob- 
vious that the localized eigenmodes are not supported by 
the structures under examination. But the calculations 
we do not present here within the framework of our sim- 
ple model exhibit the same behavior of the transmittance, 
as in Ref. 41. We have taken the geometrical parameters 
of the inclusions the same as in Ref. 41. When E is par- 
allel to Oi, the amplitude of the excited (±1,0) + SPP is 
proportional to £±i,o, which increases as the aspect ra- 
tio increases. Therefore this leads to the transmittance 
peak growth and redshift. On the contrary, when E is 
parallel to the amplitude of the excited (0, ±1)+ SPP 
is proportional to £o,±i, which decreases with an aspect 
ratio increase. This brings about the transmittance peak 
decrease and blueshift. Thus, the polarization behavior 
may be successfully described under the assumption of 
the excitation of purely interface SPPs. 

C. [1,1]+ resonance 



corresponding to the small rectangles is redshifted and 
increased as compared with the large rectangulars and 
circles as it was in the experiment. If we rotate the po- 
larization plane by ir/2 so that E is directed along g 2 , see 
Fig. 7 (b), the transmittance will depend upon the inclu- 
sion shape differently. In Fig. 9 (b) the transmittance 
is shown for different inclusion shapes both with equal 
£i for all shapes and when the modulation amplitude is 
taken from the previous case for the small rectangles. 

On the one hand, supposing the inclusion impedance, 
£j, to be constant, but at the same time changing the 
inclusion shape, we can predict the true transmittance 
through arrays of nanoparticles immersed into the film 
(each nanoparticle corresponds to an inclusion). On the 
other hand, by adjusting the inclusion impedance for dif- 
ferent shapes, we may provide the qualitative coincidence 
with experimental measurements of ELT through hole ar- 
rays. 

Consider the experiment of Ref. 41, where the effect of 
aspect ratio of rectangular holes on the transmittance of 
periodic arrays of subwavclcngth holes in optically thick 
metal films was measured. It was found that as the elec- 
tric field, E, of the incident light is directed along the 
short axes of the holes, a\, the zeroth-order transmit- 
tance peak (corresponding to [1,0]+ SPP resonance) in- 
creases and suffers a redshift when the aspect ratio of the 
holes, 02/01, is enlarged. Conversely, it was discovered 
that as E is directed along the long axes of the holes, 
d2, [1, 0]+ transmittance peak decreases and undergoes a 
blueshift when the aspect ratio of the holes is enlarged. 
The experimental measurements was also confirmed by 
strict numeric calculations. The authors attribute the 



Now examine the resonance near 1 fim (1.65 eV) in the 
vicinity of close-to-normal incidence, see Figs. 5, 6. This 
fourfold [1,1]+ resonance is consistent with excitation of 
a SB SPP at the metal-quartz interface, see Fig. 10 (b). 
A — 6 dependence is shown in Fig. 6 (square regions) for 
different polarizations and orientations of the incident 
plane. Fig. 10 (a) and (c) are the enlarged fragments of 
Fig. 6 (a) and Fig. 5. 

Now we explain the dependence of the transmittance 
upon A and 9 in the vicinity of [1, 1] + resonance, see the 
second (from above) solid square regions in Fig. 6. For p 
polarization, when E belongs to the incident plane par- 
allel to gi (ip = 0), both the pairs of the resonance waves 
corresponding to (±1,±1) and (=f1, ±1) diffraction or- 
ders are excited in the close-to-normal incidence region. 
Accordingly, in Fig. 6 (a) and in Fig. 10 (a) one can 
see the intersection of the two features that are relevant 
to two pairs of coinciding resonance curves, (1, ±1) + and 
(— 1, ±1)+, in Fig. 6 (d). Evidently, the similar dependen- 
cies take place if E is perpendicular to the incident plane 
parallel to gi (ip — 0) [s-polarization, see Fig. 6 (b)]. 
If the incident plane is oriented at ip = 7r/4, then for 
s-polarization of the incident wave [see Fig. 6 (c)] the 
magnetic field H is parallel to g 2 — gi. And the only 
excited pair of resonance waves in the vicinity of normal 
incidence corresponds to the diffraction orders (±1, =Fl) 
(resonance coefficients T_jtj Tl ). This pair provides an 
extremum of the transmittance close to the coinciding 
resonance curves (±1, ±1)+ in Fig. 6 (f). For the nondc- 
generated case, in the terms of classifying the resonance 
curves, Fig. 6 (e), there is the intersection of four features 
in the transmittance which conforms to four distinct res- 
onance curves: [1, 1] + . 
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FIG. 10: (Color online) (a) The dependence of the zeroth-order transmittance, to, upon the wavelength and incident angle 
corresponding to [1,1]+ resonance on Cav silver array of cylindrical inclusions, E||gi. In (b) the reciprocal grating is shown 
[compare with Fig. 7 (b)]. In (c) the wavelength dependence for the zeroth-order transmittance, to, and reflectance, po, (solid 
curves) along with the nonzeroth-order transmittance, ty+^o), (dotted curve) and total transmittance, rtot, (dashed curve) 
are shown. It corresponds to section of (a) along the plane, marked in grey. The inset presents the Rayleigh anomaly. The 
parameters of the array are the same as in Fig. 5. 
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FIG. 11: The squared field distribution for [1, 1]+ SPP reso- 
nance is shown at A = 1.0168/ira, 8 — 0. The parameters of 
the array are the same as in Fig. 5 



The principal difference of [1, 1]+ resonance from the 



[1,0]+ resonance considered above is that additional 
propagating waves exist in the quartz substrate. Namely, 
the waves having [1,0] diffraction orders provide the 
nonzeroth-order ELT 31 ' 32 due to the scattering of the res- 
onance waves through amplitudes £[i,oi- This scattering 
contributes additionally to the resonance width. Black 
dots, placed within the solid circle in Fig. 10 (b) indicate 
the positions for tangential components of wavevectors 
of these waves. On the other hand, the resonance [1, 1] + 
is less pronounced as compared with [1,0]+ one not only 
because of additional radiation losses. As one can see 
from Fig. 4, the amplitude which is responsible for 

resonance waves excitation within the limits of the main 
approximation, is smaller than the amplitude £[i,o], an d 
hence the SPP excitation is less effective. It should be 
noted that the SPP is excited not only owing to first- 
order scattering of the incident wave by [1, 1] resonance 
harmonics, but also due to the higher-order scattering 
processes. While the cubic processes (proportional to £ 3 ) 
can be neglected, the second-order scattering may com- 
pete with the linear one in some regions of parameter 
a/p, going beyond the approximation (19)-(22). Thus, 
for a/p > 0.35, the contribution to the amplitude of the 
resonance wave from the second order term, proportional 
to £p j , is of the the same order as the contribution from 

the linear term proportional to £[1,1] ■ 

All the propagating field harmonics contribute to the 
transmitted energy flux in the substrate media. While 
the parts of the energy flux relating to the zeroth and 
nonzeroth diffraction orders are well separated at a grate 
distance from the film, they are overlapped in the region 
sufficiently close to the film. Therefore, when making 
the experimental measurements, the detector position is 
crucial: if it is far away from the film, it fixes cither 
zeroth-order flux, tq [see Fig. 10 (c), solid curve], or 
nonzeroth one, T[ li0 ] [see Fig. 10 (c), dotted curve], de- 
pending upon detector orientation relative to the film. 
In case the detector is in the close-to-the-film vicinity, it 
fixes the total flux [see Fig. 10 (c), dashed curve]. The in- 
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terference pattern resulting from the existence of several 
homogeneous waves is quite noticeable in the far-field 
region in Fig 11 (b). The near-field structure may be 
understood in the same way as in the case of [r, 0] reso- 
nance. For instance, when E is oriented along gi+g2, the 
intensity distribution on the quartz side for strictly nor- 
mal incidence is proportional to |E+| 2 ~ sin [k( 1;1 ) t r] = 
sin 2 [(gi + g2)r], see the right-hand part of Fig. 11 (a). 
When E is oriented with respect to g x + g 2 at ip = it/ 4, 
the near-field structure is formed by the interference of 
the eigenmode for sin[(gi + g2)r] and the eigenmode for 
sin[(gi — g2)r], so that the field at the interface takes 
the form |E+| 2 ~ {sin[(gi + g 2 )r] + sin[(gi - g 2 )r]} 2 . 
This is in accordance with the left-hand part of Fig. 11 
(a). It is interesting to note that the field intensity is in- 
creased additionally in the regions where the inclusions 
are located. 



D. [1,0] resonance 

It should be emphasized, that SB SPP resonances per- 
tinent to the superstrate face provide strong nonzeroth- 
order ELT on condition that e_ < e+. It is clear from the 
geometrical reasons. Since the modulus of the wavevec- 
tors in the superstrate media are less than those of the 
substrate media, the metal-superstrate SPP wavevector 
{kgpp ~ k^feZ) corresponds to the propagating wave in 
the substrate media. This enables SPP leakage without 
scattering via the modulation. Another feature of super- 
strate resonances is far deeper minima in the zeroth-order 
reflectance as compared to the resonance on the sub- 
strate face. This is because the zeroth-order reflectance 
is strongly dependent upon the efficiency of SPP excita- 
tion. The incident light excites SPP on the superstrate 
face directly via the modulation, while for excitation on 
the substrate face, amplitude of the light decreases prior 
to the excitation caused by the tunneling through the film 
(thereby making the excitation process less effective). 

Let us discuss the SB metal-air [1,0] _ SPP resonance 
that corresponds to the transmittance (reflectance) max- 
imum (minimum) at A ~ 0.92/im in Fig. 5. The en- 
larged fragment of the wavelength-dependent zeroth- and 
nonzeroth-order transmittances are shown in Fig. 12 (a). 
While the zeroth-order transmittance profile is similar to 
that corresponding to the [1,0]+ resonance [their max- 
imal amplitudes are practically equal, being of order 
tq <~ exp(— 2$')], the reflectance has considerably deeper 
minima [cf. Fig. 12 (a) and Fig. 7 (c)], in accordance with 
the above-mentioned general property of superstrate res- 
onances. The amplitudes of nonzeroth-order outgoing 
transmitted waves exceed significantly the zeroth order 
amplitude [compare the solid curve for to and the dotted 
curve for tmq] in Fig. 12 (a)] since their amplitudes are 
of order T[i. ] ~ cxp(— 2$'). It follows from the Eq. 
(Bl). This estimation is quite universal and coincide with 
that made using the expression Eq. (24) for the simplest 
resonance. The enhanced nonzeroth-order transmission 



may be used for developing the light-splitters, since the 
waves propagate at an angle arccos(/3 r |+ v /e+~) relative to 
z axis. The unique feature of such a splitter is that all the 
four nonzeroth-order waves are linear polarized. Their 
polarization coincides with that of SPPs excited. 33 

Note that not only [1,0] diffraction orders are respon- 
sible for the nonzeroth-order ELT. As seen from Fig. 12 
(b), the diffraction orders [1, 1] correspond to the propa- 
gating waves in the quartz substrate as well. The reso- 
nance scattering contributes to these amplitudes ampli- 
tudes so that they are proportional to 6 2 , . Their inten- 
sities have been taken into account when calculating the 
total transmitted energy flux shown by the dashed line 
in Fig. 12 (a). 

It has to be pointed out that the transmittance behav- 
ior with respect to the polarization of the incident wave 
change, in the vicinity of [1,0] _ resonance, is similar to 
that in the vicinity of [1, 0]+ resonance. This is described 
in Subsection III A. (cf. the dashed square regions of A— 
and the upper solid square regions in Fig. 6.) The near- 
field distribution in X — Y plane is likewise equivalent 
to that corresponding to [1,0]+ resonance, but now the 
field is localized at the air-metal interface. This is clearly 
seen in Fig. 13 (b). Note that the far-field structure has 
the form of the interference pattern both in the air and in 
the quartz half-spaces. The pattern in the air superstrate 
is due to the interference between incident and reflected 
zeroth-order waves: we see from Fig. 12 (a) that the re- 
flectance is of order of 20% in the resonance. The pattern 
in the quartz substrate is caused by the interference be- 
tween the transmitted zeroth-order and nonzeroth-order 
waves. 

Let us illustrate the dependence of the zeroth and non- 
zeroth order transmittance upon the modulation am- 
plitude. The transmittance (reflectance) has the well- 
defined maximum at & ~ 2£/, which corresponds to 
lf[i,o]l - 0.6^^ - v /|£, see Fig. 12 (c). So, as ex- 
pected, this value of the modulation amplitude guarantee 
exactly the resonance harmonic optimal value, see discus- 
sion below Eq. (30), and, therefore, leads to the optimal 
transmittance value. 



E. Double bounday SPP 

DB SPPs for nonsymmetrically-surroundcd film ex- 
ist under a specific relation between angle of incidence, 
wavelength, structure spacing, and dielectric permittivi- 
ties of the surrounding media 32 . For instance, in the case 
of strictly normal incidence the excitation of four reso- 
nance waves in the superstrate (r = +) having multiin- 
dexes [r+ , 0] , and four resonance waves in the substrate 
(r = — ) having multiindexes [r_,0] may take place if 
K + 1 K- = r+/r_ . We do not consider this specific set of 
parameters in the present paper, but study DB SPPs in 
a symmetrically-surrounded film below. 
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FIG. 12: (Color online) In (a) the wavelength dependence for the zeroth-order transmittance, To, and reflectance, po, (solid 
curves) along with the nonzeroth-order transmittance, Tr±im> (dotted curve) and total transmittance, r to t, (dashed curve) are 
shown for [1,0]_ resonance on silver C4„ array of inclusions, E||gi, 6 = 0. In the inset the Rayleigh anomaly is shown, (b) 
The reciprocal grating [compare with Figs. 7 (b), 10 (b)]. (c) The dependence of to, and po upon the modulation amplitude, 
w = 1(0 — C0/OI- The parameters of the array are the same as in Fig. 5. 
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FIG. 13: The intensity distribution is given for [1,0]- SPP 
resonance at A = 0.9234 pm, 6 = 0. The parameters of the 
array are the same as in Fig. 5. 



IV. DB SPP RESONANCES IN THE 
SYMMETRICALLY-SURROUNDED FILM 



To the best of our knowledge, the earliest observa- 
tions of the ELT through hole arrays due to the DB 
SPP excitation were performed in Refs. 8,9. Moreover, 
in Ref. 8 the array was formed in the metal film sur- 
rounded by quartz from one face and by a liquid with 
a close-to-quartz dielectric constant from the other face, 
while Ref. 9 deals with a free-standing film. In the above 
experiments the matching of wavevectors at the opposite 
interfaces was made to excite a DB SPP. This consid- 
erably increased the ELT efficiency in zeroth diffraction 
order as compared to the SB SPP excitation for a non- 
symmetrically surrounded film. 

As known, there exist two SPP modes in the nonmod- 
ulated symmetrically-surrounded metal film, see Sec. II. 
One of them, a LR mode, corresponds to the nonsym- 
metric surface charge distribution, and higher frequency 
and Q-factor as compared to the other, SR mode. A re- 
markable feature is that the similar situation is inherent 
to the mechanical system of two coupled damping oscilla- 
tors which is nothing else but a classical analog of coupled 
SPPs in the film, see Appendix A. If the modulation is 
not too high, LR and SR modes undergo shift and widen- 
ing, but remain well-separated from one another. This 
is not exactly what is observed in the experiments using 
subwavelength hole arrays. The modulation amplitude 
of the arrays is too high to distinguish well the difference 
between LR and SR modes. 

The transmittance spectra of a symmetrically- 
surrounded film is shown in Fig. 14. One can see the four 
peaks: the first pair is 0.663 /im, 0.683 fj,m in the vicin- 
ity of [1, 1] S - L DB SPP resonance, and the second pair is 
0.918 fim, 0.929 fim in the vicinity of [1,0] 5 ' L DB SPP 
resonance. Indexes S, L comply with LR and SR modes. 
In Figs. 15 (a-c) and (e) the transmittance spectra are 
presented as functions of the wavelength and/or energy 
and the incident angle for different polarizations of the 
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FIG. 14: (Color online) Dependence of transmittance and re- 
flectance upon the wavelength for strictly normal incidence 
onto the symmetrically surrounded (by air) silver film. The 
periodical array possesses C± v symmetry, i.e. pi — pi = p — 
0.9pm, a/ p = 1/3. E is parallel to gi (to x axis, see. Fig. 1). 
The impedance of inclusions, is assumed to be £, = 2£f. 
The film thickness is equal to 2.75 skin depths. The transmit- 
tance and reflectance corresponding to the unmodulated film 
are shown by dotted lines. Vertical solid (dashed) lines indi- 
cate the wavelengths for SPPs at the unmodulated boundaries 
(Rayleigh points), one of them is shown in the inset. 



of SR mode is higher than that of LR mode, — ~ 
sinh <j>' — <j>' 

— : > 1, which leads to the higher Q-factor of 

sinh $' + $' 

LR mode as compared with SR mode. 

Second, the resonance peaks for LR and SR modes 
is strongly dependent upon the exponential index ap- 
pearing in the field approximation inside the film. In 
other words, assuming that k in Eq. (4) is complex and 
taking into account both real and imaginary parts of 
the impedance of the conductor in Eq. (4) the trans- 
mittance maxima (reflectance minima) for LR mode are 
higher (deeper) than those for SR mode. And vice versa 
if dissipation losses are neglected. As was indicated in 
the experimental paper 60 , where plasmons were stud- 
ied in the symmetrically surrounded corrugated silver 
film, "this process has a complex distance (film thick- 
ness) dependence", that is "for thin silver films (< 30 
ran), the symmetric coupled SPP is strong but extremely 
sharp, whereas the antisymmetric SPP is broad but weak. 
As the silver thickness increases, the symmetric mode 
weakens and broadens, whereas the antisymmetric mode 
sharpens and intensifies". Thus, to our opinion, the dis- 
cussed question of the resonance peaks for LR and SR 
modes needs further careful studying. 



incident wave. In Figs. 15 (d,f) the resonance curves are 
shown by solid lines. 

Consider the deviations from the normal incidence. In 
Figs. 15 (a-c) the special symmetry is shown, that is, the 
incident plane is oriented relative to gi at and 7r/4. 
This complies with coincidence of the resonance curves, 
see discussion below formula (18). It is seen in the vicin- 
ity of [1, 0] resonance shown by square regions in Fig. 15. 
In Figs. 15 (a-c) from six to two resonance features be- 
long to these regions. These features correspond from six 
to two resonance curves in Figs. 15 (d,f). Oppositely, in 
Fig. 15 (e) the transmittance contour plot is shown for 
non-specific geometry, and therefore, eight features (or, 
eight curves) belong to each vicinity according to fourfold 
DB SPP. 

Note the two features in Fig. 14. First, the widths of 
resonance peaks (dips) are less for LR modes than for SR 
modes. This is due to the higher Q-factor of LR mode 
as compared with SR mode, that can be easily explained 
in terms of field structure of SPP modes in the metal, 
^-dependence of the electric field tangential component 
for LR and SR eigenmodes is defined by Eq. (35). Since 
the amplitude of the tangential component of the elec- 
tric field is higher than the z-component within the film, 
vector E( makes the main contribution to the loss power, 
P, so that 




In view of Eq. (35), one can see that the Ohmic losses 



A. [1,0] resonance 

Now consider the [1,0] DB fourfold SPP resonance 
for strictly normal incidence. The enlarged fragment of 
Fig. 14 in the vicinity of the peaks is shown in Fig. 16 
(a). The LR and SR resonances are blueshifted as com- 
pared with the wavelengths for nonmodulated film SPPs 
marked by the vertical lines. Analogously to the SB SPP 
resonance, the amplitudes of the resonance waves are pro- 
portional to £[io], and the resonance contribution to the 
amplitude of the only propagating zeroth-order wave is 
proportional to |£[io]| 2 - 

For a SDB resonance, the transmittance and re- 
flectance have two extremums. This results from the 
excitation of LR and SR SPPs, corresponding to the ze- 
ros of the denominators Am o] in Eqs. (29), (30) (and 
Appendix B). The peak value of the transmittance 
for optimal modulation amplitude may be estimated as 
T e> ~ Co" 1 ex P( — *£''), which exceeds the transmittance 
through the unmodulated film by a factor of £^" 1 » 1. 

The distance between the two extremums of the trans- 
mittance (reflectance) may be estimated by extracting 
the difference between the wavelengths corresponding to 
SR (LR) SPPs from the dispersion relations Eq. (17). 
Using (i L j^ — e _1 / 2 \/l — (cq/^/euj L ' S ) 2 , (where we have 
replased k^vij by the SPP wavevector, q), we find for 
the frequencies uj L ' S — cq/VeB L > s , where for the LR 
mode B L = 1 - e£ tanh 2 ($/2) and for the SR mode 
B s = 1 — e£ coth 2 ($/2). In terms of wavelength we 
have for the relative difference between the wavelengths 
(X s - X L )/X ~ |£ | 2 exp(-$'), where A is of order of 



19 







-J 





Q &3 




Q 7.1 



10 



f: JO 
flfikienl a roe i.ae^remt;- 



2H 




nodc-nt angle ild^gnres) 



in -24 



nctdcnl ange (decreet!' 



-iD 



D iKJ 
ncidcntanglc Idegrccsji 



FIG. 15: (Color online) Dependence of transmittance upon the wavelength (energy) and the incident angle for different 
polarization of the incident wave and different orientations of the incident plane. The parameters of the array and the film are 
the same as in Fig. 14. (d, f) the resonance curves [see Eq. (18)]. Squares indicate the vicinity of the [1,0] resonance which 
will be considered in detail below. 



A 5 , X L . The modulation contributes additionally to this 
difference, and this contribution may be found from the 
dispersion relation An = 0, Eq.(27), see Ref 32. As 
the film thickness increases, the distance between A and 
X L decreases, see Figs. 16 (a), (c). When the resonance 
width, AA/A ~ |£ol£c>i takes the value of order of the 
splitting between the LR and SR modes, AA ~ X s — X L 
[which occurs for film thickness $' > InQ^p 1 )], the 
two maxima of the two-humped resonance curve become 
indistinguishable. 



shown in Fig. 16 (d), (e). One can see that LR (SR) 
mode is antisymmetric (symmetric) with respect to the 
tangential component of the electric field and the sur- 
face charge distribution (the surface charge sign coin- 
cides with the sign of the production n z E z , where n z 
is z-componcnt of the surface normal vector pointing out 
of the metal). 



The near-field distribution for LR and SR SPPs is 
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FIG. 16: (Color online) (a) The dependence of the zeroth-order transmittance, to, and reflectance, po, upon the wavelength 
corresponding to [1,0] resonance on the silver array of inclusions (with the symmetry Ca v ), E||gi, 9 = 0. (b) the reciprocal 
grating [compare with Figs. 7 (b), 10 (b), 12 (b)]. (c) the dependence of to upon the wavelength and film thickness for 9 — 0. (d, 
e) the distribution of tangential and normal components of the electric field at A = 0.918pm (for LR mode) and A = 0.9285/im 
(for SR mode), 9 — 0. The inclusions in (d) are shown by the dashed lines. The parameters of the array are the same as in 
Fig. 14. 



CONCLUSIONS 



We have treated analytically the resonance optical 
properties of 2D periodically-modulated optically-thick 
metal films. Explicit analytical expressions for transfor- 
mation coefficients related to any diffraction order have 
been obtained. In studying the complicated multiple res- 
onances, we have ascertained that the most of the physi- 
cal properties of ELT may be well understood in terms of 
the simplest example of a solitary resonance (when SPP is 
excited in a single diffraction order). We have examined 
and explained not only the amplitude and polarization 
dependencies of the transmittance and reflectance upon 
the parameters (angle of incidence, wavelength, tilting 
angle, film thickness, etc.), but also the field structure 
of the diffracted light. We have shown that according to 
the conception of interface SPPs excitation, the polariza- 
tion dependencies of the light diffracted by the periodical 
array may be adequately described. We have compared 
our theoretical calculations with recent experiments and 
found an excellent coincidence. We have shown that for 
hole arrays (which are in most cases used for the exper- 



imental study of ELT) it is difficult to make a distinc- 
tion between the long-range and short-rang SPP modes 
excitation caused resonance features since the radiative 
broadening exceeds the splitting between them for rather 
thick films. In this sense structures with weak modula- 
tion (say, the periodical arrays of metal nanoparticles im- 
mersed into a metal film, corrugated shallow diffraction 
gratings) are more preferable. Moreover, the efficiency of 
SPP excitation may be much higher for the structures 
with weak modulation as compared to the structures 
with holes being included. This is due to the existence 
of the optimal modulation amplitude. It was demon- 
strated that in some cases the energy flux corresponding 
to nonzeroth diffraction orders, exceeds that for zeroth- 
order. This can have a strong influence upon the results 
of experimental measurements. 



A highly instructive analogy was drawn between mo- 
tion of the monochromatic force driven system of two 
weakly coupled classical oscillators and ELT effect. 
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APPENDIX A: TWO COUPLED DAMPING 
OSCILLATORS 

The enhanced transparency and other resonant prop- 
erties of the periodically modulated metal film can be 
described by analogy with a simple mechanical system 
consisting of two weakly coupled harmonic oscillators un- 
der the action of the harmonic force applied to one of 
them see Fig. 17 (a). An oscillation amplitude of the 
second oscillator is strongly dependent on the frequency 
of the force applied to the first oscillator and can reach 
amplitudes of order of the amplitude of the first uncou- 
pled oscillator. We will concentrate on the system of two 
identical oscillators which has the analogy with the ELT 
through symmetrically-surrounded film. Lagrangian of 
this system is 



4) 



ax\X2 

+F 1 {t)x 1 + F 2 (t)x 2 , (Al) 

where dots denote time derivatives, the term ax\x 2 is 
responsible for weak coupling, a C 1, and the two last 
terms are for the forces acting on the oscillators. Adding 
to the Lagrangian —k(xf + x%)/2 we arrive at the two 
masses linked by a string with elasticity coefficient k. 
Units are such that the unperturbed frequency equals 
unity. The losses are taken into account by the dissipative 
function 63 , 



T = fii(x 1 + x 2 ) + 2/i 2 ±iX 2 • 



(A2) 



where fii > \n 2 \ > (non-negativeness of F). Below the 
dissipation is supposed to be small, \ni, 2 \ <C 1. Then 
equations of motion are 

dt dii dxi dii ' 

In normal coordinates, X^ = {x\ ± x 2 )/V2, which 
correspond to the inphasc (A_) and antiphase (X + ) os- 
cillations (symmetric and antisymmetric with respect to 
transposition 1^2 modes), we split the system (A3), 

X T + (l T a)X T + 2( f i 1 ± f i 2 )X T =F T , (A4) 

where F T = (F 1 ± F 2 )/V2. 

Solution of the corresponding homogeneous system is 



X± = A± cxp(—iil±t) , 



(A5) 



where the complex eigenfrequencies are ri=p and — O^- 
($tefL T > 0, 3?mO T < 0), 



fin 



± n 2 ) + y/l =F a - (/xi ± fi 2 ) 2 
a (^i±^ 2 ) 2 



Eigenfrequency of the inphase (symmetric) mode is lower 
than that for antiphase (antisymmetric) one, 5ief2_ < 
3?eO + if the damping is small, (/ii ± /Z2) 2 <C a. This 
is in compliance with the modes of the symmetrically- 
sandwiched film. Namely, eigenfrequency of the LR (an- 
tisymmetric) SPP is lower than that for SR (symmetric) 
one, see discussions in the Subsection IV A. In turn, the 
corresponding decrement, T T ~ /ii± / u 2 , is higher (lower) 
for /i 2 > (p 2 < 0). Note that the case with /j, 2 > holds 
for different active losses for symmetric and antisymmet- 
ric polariton modes. This difference is due to the fact 
that the high-frequency mode has a lower mean value of 
the electric field in the film, and, hence, lower losses and 
a higher Q-factor. 

The specific solution of the dynamic equations (A4) 
for harmonic driving forces being applied, F\ <2 = 
/i,2 exp(— iflt), is 



± ! 



(A7) 



X±(t) = A±(Q)exp(-int), A± 
where f± = (fi =F /2)/V2, 

r»± = x»±(0) = -(n-n±)(n + n*±). (A8) 

Thus, displacements of the forced oscillations are Xi(t) = 
a,i(Q) exp(-iQt), i = 1,2, 



ai, 2 (0) = 



~7T~ 



h + h ± h 

2V_ 



2V A 



(A9) 



~ 1 =F 



i(Mi±M2)- (A6) 



Let excitation of the system be caused by a driving 
force applied to the first oscillator, fi = 4, f 2 = 0. This 
is analogous to the diffraction problem, where only one 
wave is incident onto a face of the metallic film. If the 
dissipation is small, (^1 ±^ 2 ) 2 -C a, then there exist two 
pronounced maximal magnitudes of a\ y2 in the vicinity of 
the frequencies Q = where ai t2 (n'_) — 2/2?_(f2'_) ~ 
i/T-, a h2 (fl' + ) ~ ±2/D+(0' + ) ^'±i/T+. The splitting 
of the resonance maxima is approximately a, and the 
up- (low-) frequency resonance widths are SCl± = 2T±. 
As a decreases, the distance between resonance peaks 
diminishes, and for a <~ T + +r_ they overlap, see Fig. 17 
(b) , (c) . The mechanical coupling factor a is analogous to 
that arising in the EM problems: for the splitting of the 
SPP modes in the film this factor is exp(— <E>'), compare 
Fig. 17 (c) and Fig. 16 (c). 

For /Hi = fj, 2 = the amplitude of the enforced oscilla- 
tor vanishes at the unperturbed eigenfrequency, = 1. 
This well-known effect of the resonance damping (which 
is extensively used in shipbuilding) also shows an analogy 
with the ELT problem, where Fano minima arise in the 
amplitudes of the fields. 

In a similar way one can draw an analogy for the 
nonsymmetrically-sandwiched film. In that case the me- 
chanical system consist of two oscillators of different 
masses. 

Thus, some physical properties of the ELT phenomena 
are quite general. At least, the similarities with the sim- 
plest classical mechanical systems may be established. 
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FIG. 17: (Color online) (a) The geometry of the problem and the mechanical analogy, (b) The dependence of "free" oscillator 
amplitude upon the period, 1/fi. The coupling parameter and the damping are: a — 0.05, 0.03, 0.015, 0.01 for the solid, 
dash-dotted, dotted, dashed curve respectively; Fi = 0.004, T2 = 0.001 for all curves, (c) Dependence of the "free" oscillator 
amplitude on the period, 1/fi, and coupling parameter a. 



APPENDIX B: TRANSFORMATION 
COEFFICIENTS FOR FOURFOLD RESONANCE 
IN CASES OF C 2v AND C iv SYMMETRIES 

In the case of Ci v symmetry which implies the follow- 
ing properties of modulation impedance: 

£ni,n2 ni,n2 Cni,— ni £— m,— n2 3 

the resonance matrix possesses the symmetry properties. 
The resonance TCs in the vicinity of normal incidence for 
[r, 0], or [r, 0] T resonance are related as T^ 17 = —T_~q, 

Tq t a — —Tq _°. This leads to the essential simplifica- 
tion of Eqs. (19), (22): the resonance matrix 8x8 can 
be reduced to a matrix 2x2. Thus, the TCs are written 
as 



Tj} aa '= 5 M ,o5^'T T F W + ^ k (cosM>)-^ 



n 



[h\n + T \n - ^(cosh^" 1 ] , (B3) 



where AZj-, 



l-a 1 + " 1 + ° 

2ct't— e r 4 e_ 4 



P-\o PtW- If the ad- 
ulation has C^ v symmetry, the zeroth-order TCs allow 
additional simplifications 



T. 



,t|— a 



n 



r— 1 



0t\K ~ T TC (cosh$ 



where 1Z = (r, 0), or 1Z = (0, r), 



(cosh$) * 



Ar = P+\nP-\n - T^cosh 2 $, 
0r\n = /3 r | TC tanh$ + Co + 6rc + Gji, 
= £o + + + G n , 



AT 



-n 



(Bl) 



(B2) 



T, 



o 



x /3t|k4 
x (cosh$) 



l + T 



T [r>0] )(cosh$) T - 1 



t,i i<tct = o. 



(B4) 



Notice that the zeroth-order TCs are diagonal in po- 
larization, i.e. the polarization of zeroth-order re- 
flected/transmitted wave is the same as that of the inci- 
dent wave. 
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For a modulated film there can exist coupling between SB 
SPPs caused by diffraction. Such DB dressed SPPs cor- 
respond to specific relations between permittivities of the 
dielectric media. Formally they correspond to simultane- 
ous vanishing of two determinants, ||b^j ' CT || and ||6^/' CT ||, 



with M' 7^ M, and in calculating the matrix determinant 

-^mm* ( m ^ ne lowest order in the modulation) one has to 
take into consideration the products of the determinants of 
the diagonal in diffraction orders submatrices b^}," along 
with the terms corresponding to the submatrices d^^, CT , 

jTT \oO 

a M'M ■ 

In case of the periodicity formed by the surface relief cor- 
rugation, the situation changes principally: coupling does 
exist even when SPPs propagate perpendicular to one an- 
other. 



